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ABSTRACT 

We monitored five high-luminosity quasars with ^ > 10 45 er S s 1 a t 0.4 < z < 0.6 to 

measure flux variability of the Mgll A2798 line and explore feasibility of reverberation mapping using 
Mgll. Over the two year monitoring program, imaging data were obtained with the A Noble Double- 
Imaging Camera on the 1.3-m telescope at the Cerro Tololo Inter- American Observatory (CTIO), 
while spectroscopic data were obtained at the same night with the R-C spectrograph on the 1.5-m 
telescope at the CTIO. By performing differential photometry using available field stars in each quasar 
image, we measured variability - 10%-24% peak-to-peak changes and 3%-8% rms variations - in the 
B band, which includes flux changes in the rest-frame UV continuum (^2500A- ~3600A) as well as 
the Mgll line. Utilizing photometric measurements for spectroscopic flux calibration, we measured 
the Mgll line flux and the continuum flux at 3000A from each single-epoch spectrum. Four objects 
showed Mgll line flux variability with 23%-50% peak-to-peak changes and 8%-17% rms variations 
over 1-1.5 year rest-frame time scales, while one object showed no Mgll flux variability within the 
measurement error (<5%). We also detected 4%-15% rms variations of the Mgll line width for all 
five objects. With synchronous observations for photometry and spectroscopy, we demonstrated the 
feasibility of the Mgll line reverberation mapping for high- luminosity quasars at intermediate redshift. 
Subject headings: galaxies: active — galaxies: nuclei — quasars: general 



1. INTRODUCTION 

Determining black hole mass is a fundamental step to 
understand the physics of the AGN accretion disk (e.g. 
Bonning et al. 2007; Davis et al. 2007) as well as cosmic 
evolution of black holes and galaxies (e.g. Di Matteo et 
al. 2005; Croton et al. 2006; Woo et al. 2006; Treu et al. 
2007). Active galaxies with broad emission lines provide 
a way of estimating black hole masses, assuming that the 
motion of the broad-line regions (BLR) is dominated by 
the gravitational potential of the central source. Thus, 
virial mass can be determined from the dynamics of the 
BLR. The width of broad-emission lines are used in es- 
timating the velocity scale, while the size (Rblr) can 
be determined by measuring the time lag between the 
continuum flux changes (presumably originating from an 
accretion flow) and the corresponding flux changes in the 
broad-emission lines (Blanford & McKee 1982). 

So far, reverberation masses have been determined for 
about three dozens of local Seyfert galaxies and PG 
quasars, which are limited to redshift z < 0.4 and op- 
tical luminosity L 5100 ^ < 10 46 erg s _1 (Wandel et al. 
1999; Kaspi et al. 2000, Peterson et al. 2004), while 
there are numerous AGNs with higher optical luminosity 
at z > 0.4. Thus, it is crucial to increase the sample size, 
especially including higher luminosity and higher redshift 
AGNs, in order to validate the virial technique and un- 
derstand its uncertainties. The size-luminosity relation 
- correlation between the measured time lag (BLR size) 
and the continuum luminosity (Kaspi et al. 2005; Bcntz 
et al. 2006) - provides a simple recipe for black hole mass 
estimation without reverberation mapping (e.g. Woo & 
Urry 2002; Vestergaard 2002; Netzer & Trakhtenbrot 
2007; McGill et al. 2007). Therefore, measuring time 
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lags for higher luminosity AGNs is also necessary to val- 
idate the size-luminosity relation and avoid extrapola- 
tions in estimating black hole masses for high luminosity 
and high redshift quasars (e.g. McLure & Dunlop 2004; 
Vestergaard 2004). 

Depending on the redshift of the source, different broad 
emission lines can be used for reverberation studies. For 
local AGNs, the best studied Balmer lines, such as H/3 
and Ha, showed clear correspondence to the continuum 
changes (e.g. Peterson et al. 2004). At z > 2, CIV 
A1549 and Lycv lines are redshifted in the optical range 
and become useful for reverberation studies with optical 
spectrographs. Recently, Kaspi et al. (2007) showed fea- 
sibility of reverberation mapping for quasars at z > 2, 
reporting CIV line flux variability for six quasars and a 
tentative time lag measurement for one object. In the 
same study, however, the Lycv line showed no flux vari- 
ability within 5% measurement error. In the local uni- 
verse, time lag measurements with CIV line have been 
reported for several Seyfert galaxies (see Peterson et al. 
2004, for a summary). 

At intermediate redshift, 0.7 < z < 2, the Mgll A2798 
line becomes virtually the only available broad line for 
estimating virial black hole masses. However, there has 
been only a handful of studies on Mgll line flux variabil- 
ity, and it is not clear whether the Mgll line shows the 
same reverberation as other broad lines. For example, 
Clavel et al. (1991) reported that Mgll line flux variabil- 
ity of NGC 5548 was much lower than flux changes in 
other broad lines and did not show clear correspondence 
to UV continuum flux changes. In the case of NGC 3516, 
the Mgll line flux was almost constant in 1996 although 
other lines and the continuum showed large flux varia- 
tions (Goad et al. 1999). In contrast, time-lag measure- 
ments were consistent between Mgll and Balmer lines for 
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NGC 3783 (Reichcrt et al. 1994) and NGC 4151 (Met- 
zroth ct al. 2006). At high redshift, 1 < z < 2, Trevese 
et al. (2007) studied two high-luminosity quasars and de- 
tected Mgll line flux variability for one object. Clearly 
flux variability of the Mgll line has to be better con- 
strained to test feasibility of reverberation mapping, es- 
pecially for intermediate-to-high redshift quasars. 

We initiated a monitoring program to study Mgll line 
variability and explore feasibility of reverberation map- 
ping for bright quasars with the 1.5-m and 1.3-m tele- 
scopes at the Cerro Tololo Inter-American Observatory 
(CTIO) when these telescopes became available in ser- 
vice mode through the Small and Moderate Aperture 
Research Telescope System (SMARTS) consortium 1 in 
2003. In this paper we report our imaging and spec- 
troscopic study on the Mgll line variability of five high- 
luminosity quasars at 0.4 < z < 0.6. The goal of the 
paper is to investigate variability of the Mgll line flux 
and width. The paper is organized as follows. In § 2 
we describe observations and data reduction. In § 3 we 
describe our measurements of B band luminosity vari- 
ability and Mgll line flux variability. In § 4 we discuss 
our results and their implications. A standard cosmol- 
ogy is assumed where necessary (H =70 kms _1 Mpc _1 , 
O m = 0.3 and Q A = 0.7). 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Sample Selection and Experimental Design 

To study feasibility of reverberation mapping with the 
Mgll line, a sample of bright southern quasars was se- 
lected from the Calan- Tololo Survey catalogue (Maza et 
al. 1993, 1995, 1996) with a redshift range, 0.4 < z < 1, 
where no reverberation study has been reported. Twelve 
bright objects were initially selected based on the ob- 
servability of more than 6 months per year. Here we 
present five best observed objects with high quality spec- 
tra and strong Mgll lines, observed at more than 5 dif- 
ferent epochs for spectroscopy and photometry Target 
properties of these objects are listed in Table 1. The 
other seven objects were excluded in this study due to 
the fact that observed spectra have lower signal-to-noise 
ratio (< 10), not enough to measure line flux variations 
with ~5% measurement errors, and/or the number of ob- 
served epochs with both photometry and spectroscopy is 
typically less than 5. 

We designed spectroscopic and photometric observa- 
tions to be carried out at the same night so that spectro- 
scopic flux calibration can be accurately obtained using 
photometric measurements, assuming that there is no 
intra- night flux variations (see §3.2). However, because 
of the difficulty of scheduling, there are several epochs 
when spectroscopy and photometry were carried out at 
different nights. We excluded data from those epochs for 
the analysis of the Mgll line flux variations (see Table 
2). 

2.2. Photometry 

Photometric and spectroscopic observations were per- 
formed in service mode, respectively, with the CTIO 1.3- 
m and 1.5-m telescopes operated by the SMARTS con- 
sortium between March 2003 and February 2005. Due to 

1 http://www.astro.yale.edu/smarts 
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Fig. 1.— B band (rest-frame ~2500A-~3600A) light curves 
of all five objects. Pcak-to-peak variations are 9%-24% and rms 
variations are 3%-8%, calculated in logarithmic scale. 

the observability and limited telescope time, ~10 single- 
epoch data were obtained for given objects during the 
monitoring program. Table 2 shows journal of observa- 
tions. 

For photometry, we used the A Noble Double- Imaging 
Camera (ANDICAM), a dual-channel optical IR imager 
at the CTIO 1.3-m telescope. The optical CCD with an 
excellent cosmetics provided ~6x6 arcminute 2 field of 
view and 0.369" per pixel resolution. Each target quasar 
was roughly centered in the field of view to cover as 
many as possible field stars for differential photometry. 
We used the broad B band filter to match the observed 
wavelength range in spectroscopy, in order to do spectro- 
scopic flux calibration using photometric measurements. 
The typical exposure time was 10 minute, yielding signal- 
to-noise ratio (S/N) > 100 on the target quasars. Bias 
and dome flats were taken in the afternoon before each 
night. A set of standard stars were also taken during 
the night. Initial data reduction, including bias subtrac- 
tion and flat-fielding, has been automatically performed 
through a pipeline by the Yale SMARTS team. 

2.3. Spectroscopy 

We used the R-C spectrograph with the 600 line mm -1 
grating centered at ~ 4500A. The observed spectral 
range was ^3500-5300A, covering the Mgll line (rest 
wavelength ~ 2798A) near the center in the observed 
spectral range. The pixel scale corresponded to 1.47 
A and the spatial scale along the slit was l'/3 per pixel. 
We chose the 2" wide slit to include most of the AGN 
light within the slit under the typical sub-arsecond see- 
ing conditions at the CTIO. The spectral resolution mea- 
sured from arc lines was ~ 3.8A (FWHM), yielding Gaus- 
sian velocity of 110 km s _1 at the center of the spectral 
range. Total exposure time for each object was 3600 sec- 
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TABLE 1 
Target properties 



(1) Name 


(2)z 


(3) RA (J2000) 


(4) DEC (J2000) 


(5) B 


(6) log M BH 


(7) L/Lfidd 


CTS300 


0.352 


11 28 18.62 


-13 19 28.7 


17.2 


8.74 


0.12 


CTS306 


0.456 


11 59 41.02 


-19 59 24.7 


17.2 


8.87 


0.13 


CTS536 


0.500 


14 01 49.27 


-24 55 30.3 


16.5 


8.81 


0.40 


CTS582 


0.483 


22 40 18.84 


-52 31 52.4 


17.8 


8.75 


0.13 


CTS667 


0.557 


11 39 10.70 


-13 50 43.6 


16.7 


8.62 


1.12 



Note. — Col. (1): Target ID. Col. (2): Redshift measured from the Mgll line. Col. (3): RA. Col. (4): DEC. Col. (5): mean B 
magnitude. Col. (6): black hole mass estimated using the Mgll FWHM and the continuum luminosity at 3000A from the mean spectra. 
Col. (7): Eddington ratio calculated using black hole mass and 5.9x L • as bolometric luminosity. 





TABLE 2 






Journal of observations 




Date 


Photomcty 


Spectroscopy 




Targets 


Targets 


(1) 


(2) 


(3) 


2003 Mar 12 


CTS306, CTS536 


CTS306, CTS536 


2003 Apr 4 




CTS306 , CTS536 a 


2003 Apr 5 


CTS306, CTS536 




2003 May 7 


CTS306, CTS536 


CTS306, CTS536 


2003 Jun 17 


CTS582 


CTS536 a , CTS582 


2003 Jun 18 


CTS306, CTS536 


2003 Jun 19 


CTS306, CTS536 




2003 Jun 20 


CTS306, CTS536, CTS582 




2003 Jun 25 


CTS306 


CTS306 


onnQ Tul on 
zUUo Jul ZU 


PTCQHfi r^HTC^Qfi PTCc;CO 


OlooUD, OlaOoD 


zuuo ocp zo 




2003 Oct 21 


CTS582 


CTS582 


2003 Nov 12 


CTS582 


CTS582 


2003 Dec 16 


CTS582 


CTS582 


2003 Dec 19 


CTS300 




2003 Dec 20 


CTS667 




2003 Dec 29 




CTSSOO 6 


2003 Dec 30 




CTsee? 6 


2004 Jan 20 


CTS300, CTS667 


CTS667 


2004 Jan 21 


CTS306, CTS536 


CTS306, CTS536 


2004 Feb 23 


CTS536 


CTS536 


2004 Mar 2 


CTS300 


CTS300 


2004 Mar 16 


CTS667 


CTS667 


2004 Mar 31 


CTS300, CTS536 


CTS300, CTS536 


2004 Apr 15 


CTS300, CTS306, CTS536, CTS667 


CTS300, CTS306, CTS536, CTS667 


2004 Apr 16 


CTS582 


CTS582 


2004 May 13 


CTS582 




2004 Jun 22 


CTS536 


CTS536, CTS582 i > 


2004 Jul 14 


CTS536 


CTS536 


2004 Aug 11 


CTS582 




2004 Sep 9 


CTS582 




2004 Nov 16 




CTS582 6 


2004 Dec 18 


CTS582 


CTS582 


2004 Dec 20 


CTS667 


CTS667 


2004 Dec 21 


CTS306 


CTS306 


2004 Dec 22 




CTS300 6 


2005 Jan 12 


CTS667 


CTS667 


2005 Jan 13 


CTS306 




2005 Jan 14 


CTS300 


CTS300 


2005 Jan 15 




CTSSOe 6 


2005 Feb 8 




CTS300 6 


2005 Feb 9 




CTSee? 6 



References. — a) these spectra were flux calibrated with photometric measurements from the previous or following night, b) these 
spectra were excluded in the analysis of Mgll line flux variability due to the lack of photometric data 
Note. — Col. (1): Observing date. Col. (2): Photometry targets Col. (3): Spectroscopy targets 
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ond, divided into 3 exposures in order to remove cosmic 
rays, yielding S/N ~ 10-30 per pixel on the continuum. 
Internal lamp flats were taken in the afternoon to cor- 
rect pixel-to-pixel variations. A set of spectrophotomet- 
ry standards were observed during each night for initial 
flux calibrations (see §3.2). 

We performed the standard data reduction including 
bias subtraction, flat fielding, wavelength calibration, 
spectral extraction, and flux calibration using a series 
of IRAF scripts. He and Ar arc lines were used for 
wavelength calibration. One-dimensional spectra were 
extracted for maximal S /N with a typical extraction ra- 
dius of 3 pixels, corresponding to ~ 3.9". Individual 
spectra were combined with inverse-variance weighting 
to create the final spectrum for each single-epoch. For 
each night, we obtained sensitivity function using spec- 
trophotometric stars. The flux scale is consistent within 
1% rms variation over the monitoring period. 

3. MEASUREMENTS 

In this section, we present broad-band flux variability 
measured from differential photometry (§3.1). Then we 
show our spectroscopic analysis for Mgll line variability. 
First, we demonstrate spectroscopic flux calibration with 
photometric measurements (§3.2). In §3.3 we present 
how we subtract the continuum and measure the Mgll 
line flux. We report Mgll line flux variability in §3.4, 
and Mgll line width variability in §3.5. 

3.1. B band flux variability from photometry 

To measure the broad band flux variability, we per- 
formed differential photometry using stars in the field of 
each quasar, as demonstrated in the paper by Woo et al. 
(2007). More than 10 stars were initially used to check 
variability and photometric errors. Aperture photome- 
try of comparison stars and quasars was performed using 
the PHOT task in IRAF, with a 10 pixel (corresponding 
to 13") aperture and with an annulus between 15 and 20 
pixels for sky subtraction, based on our initial tests with 
images of various quality seeing. Photometric errors of 
individual stars were typically less than 1%. We chose 
three to five unsaturated stars brighter than each quasar 
for differential photometry. First, we calculated the mag- 
nitude difference between each comparison star and the 
target quasar. Second, we normalized the magnitude dif- 
ference to zero by subtracting the mean difference over 
all observed epochs. Photometric errors of comparison 
stars were added in quadrature and then the photomet- 
ric error of the target quasar was added on top of it, 
yielding total errors of ~1% for relative flux changes. 

In Figure 1 we show relative B band flux variability. 
All five quasars showed significant flux variability in the 
B band. We find peak-to-peak variations in the range of 
9%-24% and rms variations 2 in the range of 3%-8%. We 
note that since the Mgll line flux is also included in the 
B band photometry, B band flux variability cannot be 
considered as the continuum (such as 3000A) variability 
unless the Mgll line flux is constant (see §3.4 for more 
details). Nevertheless, the detected B band flux varia- 
tions clearly demonstrate that our sample quasars have 
UV continuum variability. 

2 We calculate rms variations using fluxes in logarithmic scale 
throughout the paper. 
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Fig. 2. — A series of continuum subtracted spectra around the 
Mgll line of CTS300. The Mgll lines are modeled with a multi- 
Gaussian profile (red). Observing dates (JD-2450000) are denoted 
in the upper right Conner of each panel. 

3.2. flux calibration using photometric measurements 

The flux uncertainty of the ground-based spectroscopy 
is typically larger than 10% due to the slit losses, sky 
transparency, and various seeing conditions. In rever- 
beration mapping studies, non-varying narrow emission 
lines, such as [O III] originated from far extended low- 
density regions, are used as an internal flux calibrator. 
Thus, flux variations of broad lines can be measured by 
normalizing each spectrum to a constant flux of the nar- 
row line (e.g. Peterson et al. 1998; Woo et al. 2007). In 
the vicinity of the Mgll line, there is no strong narrow 
line useful as an internal flux calibrator. The [Oil] line 
at 3727A is far from Mgll in spectral range and much 
weaker than Mgll, not desirable as an internal calibra- 
tor. However, since the Mgll line is much stronger than 
H/3, measuring its relative variation can be robustly done 
even without an internal flux calibrator as we demon- 
strate below. 

Instead of relying on the precision of flux calibration 
with spectrophotometric stars, we normalize each single- 
epoch spectrum to match photometric measurements ob- 
tained at the same night. Using B band response func- 
tion, we calculate synthetic magnitude from the observed 
spectrum and then normalize the spectrum by multiply- 
ing a scale factor in order to match the B band mag- 
nitude measured from photometry. The amount of flux 
calibration with a scale factor was typically 10%-20% of 
the original flux. 

We note that spectroscopy and photometry observa- 
tions were not performed at the exact same time. Thus, 
spectroscopic flux calibration is uncertain within the level 
of intra-night variations, which seem to be negligible for 
Seyfert 1 galaxies and quasars based on the study by 
Webb & Malkan (2000). Due to the limited telescope 
time for photometry, we have not attempted to measure 
intra-night variations for individual objects, which can 
be tested in the future. 

3.3. Mgll line and continuum flux measurements 
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Using the flux calibrated spectra, we measure the flux 
and width of the Mgll line in several steps as similarly 
done in other studies (e.g. McGill et al. 2007). First, 
we fit the continuum under Mgll with a straight line us- 
ing the average flux around ~2700A and ~2900A with 
40A wide windows. We linearly interpolate these two 
continuum points to define an underlying continuum. 
Second, we measure the Mgll line flux by integrating 
the flux using the continuum subtracted spectrum. We 
note that defining a true continuum under the Mgll line 
is not straight forward because of the presence of broad 
Fell emission blends (small blue bump). Hence, depend- 
ing on how the underlying continuum is defined, the mea- 
sured Mgll line flux can vary. Nevertheless, as far as we 
keep the same windows for determining the continuum, 
we can measure the relative changes in the Mgll line 
flux between two different epochs. We did not subtract 
the broad Fell emission, which could contribute small 
amount of flux to the wings of the Mgll line. However, 
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the relative flux changes in Mgll will not be significantly 
affected by the contribution of the Fell emission. 

The random error on the Mgll line flux measurements 
is typically 3% level (Table 3). By adding ~2% system- 
atic error of the flux calibration to the random error, 
we take 5% as an upper limit to the total error of the 
Mgll line flux measurements. A series of continuum- 
subtracted spectra for individual objects are presented 
in Figures 2-6, where variation of the Mgll line flux can 
be qualitatively detected. 

From each spectrum the monochromatic continuum 
flux was measured at 3000A by averaging the flux within 
40A window. We note that the size of the window does 
not affect the continuum variability. Typical errors on 
the continuum flux measurements were 2%-3% depend- 
ing on the S/N of each spectrum (see Table 3). 

3.4. Flux variability of the Mgll line and continuum 

In Figure 7 we present light curves of the Mgll line flux 
(left) and the continuum flux at 3000A (right). All five 
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Fig. 7. — Variability of the Mgll line flux (left panels) and the 
continuum flux at 3000A (right panels) for the sample. 

quasars showed continuum (3000A) flux variability, with 
peak-to-peak variations of 17%-38% and rms variations 
of 6%-15%. Continuum flux variation shows consistent 
pattern as B band flux variation in the case of CTS536, 
for which the Mgll line flux has been almost constant 
during the monitoring period. However, for the other 
objects, the continuum flux variations are different from 
B band flux variations depending on the relative contri- 
bution of the Mgll line flux. 

In the case of Mgll line flux, we detected 23%-50% 
peak-to-peak changes and 8%-17% rms variations for 
four objects. In contrast, one object, CTS536, did not 
show clear variation within the measurement error 5%. 
The Mgll line flux variation of these four quasars is sim- 
ilar to that of local Seyfert galaxies and high redshift 
quasars (e.g. Clavel et al. 1991; Trevese et al. 2007), 
implying that intermediate redshift quasars have similar 
characteristics of the Mgll line flux variability, although 
the actual amplitude is highly variable depending on in- 
dividual objects and the monitoring time base-line. 

3.5. Line width variability 

To measure the Mgll line width, we fit the line profile 
with multi-Gaussian models and determine the best fit 
model based on x 2 statistics. Fitting the line profile is 
necessary for low quality spectra, for which FWHM can 
not be easily defined on the data. We measured FWHM 
on the best fit model as demonstrated in Figures 2-6 
(red lines). In the case of CTS667, Fell contamination 
is present in the wings of the Mgll line profile and our 
multi-Gaussian model also fits the broad Fell component. 
However, FWHM of the Mgll line is not significantly af- 
fected by the shape of the wings. Thus, we simply ig- 
nored the Fell component although this will increase the 
uncertainties of Mgll line width measurements, which is 
dominated by systematic errors from continuum subtrac- 
tion and assumed to be ~5%. 

In Figure 8 we compared measured FWHM with the 
continuum luminosity for each object. The rms variation 
of the Mgll line width, ranging from 4% to 15% (see Ta- 
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Fig. 8. — Mgll line width variability. 5% error was assumed for 
the FWHM measurements. 

ble 3), is similar to the rms variation of continuum flux, 
suggesting a photoionization origin of the Mgll line. The 
reason why we do not clearly see a correlation between 
the line width and continuum luminosity, as expected 

1/2 

from a photoionization model (e.g. FWHM oc L ' „), 

3000A 

is probably due to the combined effects of measurement 
and systematic errors on the line width and relatively 
small flux variations of the continuum. 

3.6. Black hole mass estimates 

We estimated black hole mass using the virial tech- 
nique based on the empirical size-luminosity relation of 
the reverberation sample (Kaspi et al. 2005; Bentz et al. 
2006). Combining the Mgll line FWHM and the con- 
tinuum luminosity at 3000A, measured from the mean 
spectra of each object, black hole mass was determined 
with the empirical formula from McGill et al. (2007): 



M BH = 6.77 



( L 



3000 



V io 44 



ergs" 



°' 47 / FWHM Mg ii N 1 
V1000 km s" 1 



Al 



0- 
(1) 

The estimated black hole mass of the sample quasars 
ranges between 10 8 6 and 10 8 9 M (see Table 1). By 
multiplying a factor of 5.9 to the continuum luminos- 
ity at 3000A (McLure & Dunlop 2004), we estimated 
bolometric luminosity and computed the Eddington ra- 
tio. The determined Eddington ratio ranges from 0.1 to 
1.1 as tabulated in Table 1, indicating our sample quasars 
are highly accreting luminous quasars. 

4. DISCUSSION AND CONCLUSIONS 

We monitored and measured variability of 5 high- 
luminosity quasars at 0.4 < z < 0.6. From the pho- 
tometry, we detected 10%-24% peak-to-peak flux vari- 
ations and 3%-8% rms variations in the B band which 
corresponds to the rest-frame 2500A-3600A range. Af- 
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TABLE 3 
Variability measurements 



Target 


N 


N 


Mgll 


rms 


error 


L3000 


rms 


error 


F 


rms 


F 


L3000 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


CTS300 


8 


6 


.17 (50%) 


.067 (17%) 


.013 (3%) 


.08 (20%) 


.028 (7%) 


.011 (3%) 


5051 


.035 (8%) 


5197 


45.15 


CTS306 


10 


9 


.16 (45%) 


.054 (13%) 


.013 (3%) 


.12 (32%) 


.032 (8%) 


.011 (3%) 


5494 


.035 (8%) 


5540 


45.30 


CTS536 


11 


11 


.04 (10%) 


.012 (3%) 


.011 (3%) 


.07 (17%) 


.027 (6%) 


.008 (2%) 


4020 


.033 (8%) 


4101 


45.74 


CTS582 


10 


6 


.13 (35%) 


.034 (8%) 


.017 (4%) 


.09 (23%) 


.030 (7%) 


.013 (3%) 


5102 


.063 (15%) 


5144 


45.19 


CTS667 


7 


6 


.09 (23%) 


.036 (9%) 


.012 (3%) 


.14 (38%) 


.060 (15%) 


.008 (2%) 


2907 


.015 (4%) 


2862 


45.99 



Note. — Col. (1): Target ID. Col. (2): Number of observed epochs for spectroscopy. Col. (3): Number of epochs with flux calibration 

using available photometry. Col. (4): peak to peak variation of Mgll line flux in logarithmic scale. Col. (5): rms variation of Mgll line 

flux in logarithmic scale. Col. (6): mean error on Mgll line flux measurement in logarithmic scale. Col. (7): peak to peak variation of 

continuum luminosity at 3000A in logarithmic scale. Col. (8): rms variation of continuum luminosity at 3000A in logarithmic scale. Col. 

(9): mean error on continuum luminosity measurement in logarithmic scale. Col. (10): mean FWHM of the Mgll line in km s~ 1 . Col. 

(11): rms variation of FWHM in logarithmic scale. Col. (12): FWHM of the Mgll line in km s _1 , measured from mean spectra. Col. (13): 

L„„„„ °. m er f? s_1 > measured from mean spectra. 
3000A ° ' 1 

ter spectroscopic flux calibration using photometric mea- 
surements, we measured the Mgll line flux and the con- 
tinuum flux at 3000A from the observed spectra. Mgll 
flux variations were detected for 4 objects with 23%-50% 
peak-to-peak variations and 8%-17% rms variations. One 
object did not show clear variations in the Mgll line flux 
within the measurement error (~5%) over the two year 
monitoring period. By measuring FWHM of the Mgll 
line, we measured variability of line width for the same 
quasars, which showed 4%-15% rms variations. 

These variability detections are consistent with a pic- 
ture that the Mgll line is variable as other broad 
lines, such as Balmer lines or high ionization lines (Lya 
and CIV). The Mgll line of these intermediate redshift 
quasars probably originates in the same region as H/3, as 
shown in the cases of NGC 3783 and NGC 4151, for which 
similar time-lags were measured from H/3 and Mgll lines 
(Reichert et al. 1994; Peterson et al. 2004; Metzroth ct 
al. 2006). It is possible that for some particular objects 
Mgll line emitting region is more extended than Hf3 line 
emitting region, and hence Mgll line variability may not 
be detected in short time scales. This scenario was sug- 
gested for NGC 3516, which showed no Mgll line flux 
variations (Goad et al. 1999). However, it is not clear 
whether NGC 3516 is an abnormal object with an ex- 
tended Mgll line region, or substantial fraction of quasars 
have a similarly extended Mgll line region. Monitoring 
flux variations of the Mgll line together with the H/3 line 
using a larger sample of AGNs is necessary to properly 
address these issues. 

Although we did not attempt to do reverberation anal- 
ysis due to the short time baseline with insufficient num- 
bers of epochs, it is shown that the Mgll line can be used 



for reverberation mapping for even luminous quasars at 
intermediate redshifts. Compared to the local Seyfert 
galaxies, a longer-term monitoring program is necessary 
for these objects due to the 1+z effect on the observed 
time scale as well as larger intrinsic reverberation time 
scales for more luminous quasars. The estimated time- 
lags for the sample quasars, based on the continuum 
luminosity at 3000A and the size-luminosity relation 
(McLure & Dunlop 2004), range from 74 to 217 light- 
days. Thus, a few year monitoring program with syn- 
chronous photometry and spectroscopy can reveal the 
reverberation time scales of the Mgll line for the sample. 
Small telescopes (l-3m) can be very useful tools for Mgll 
line reverberation mapping since Mgll is much stronger 
than H/3 which has been generally used in reverberation 
mapping of local AGNs. Regular monitoring and service 
observation is a key for a success in future reverberation 
mass measurements. 
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